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APPLICATION  OF  THE  METHOD  OF  HARMONIC  LINEARIZATION  FOR  INVESTIGATING 
SLIP  MODES  IN  A  DIGITAL  TRACKING  SYSTEM  WITH  VARIABLE  STRUCTURE 

Ye.Ya.  Ivanchenko,  V.A.  Kamenev 

At  the  present  there  are  not  any  accurate  methods  for  solving 
a  nonlinear  differential  equation  of  any  order.  Therefore,  used  very 
frequently  are  approximation  methods  to  study  nonlinear  control  sys¬ 
tems,  which  allow  obtaining  results  with  an  accuracy  sufficient  for 
practice.  In  the  proposed  work  the  problem  is  stated  concerning  the 
application  of  the  widely  known  method  of  harmonic  linearization  of 
Ye.P.  Popov  for  investigating  slip  methods  in  a  system  with  variable 
structure  (SPS). 

A  very  important  characteristic  of  the  SPS's  is  the  fact  that 
they  allow  ensuring  a  high-quality  control  of  the  free  movement  of 
objects  with  constant  parameters,  and  in  the  case  of  a  change  in  the 
parameters  of  an  object  over  a  wide  range  they  allow  constructing  a 
system  which  is  barely  sensitive  to  these  changes. 

In  reference  [1]  a  detailed  description  is  given  of  the  principle 
of  construction  of  the  SPS,  and  methods  of  the  phase  plane  are  widely 
used.  Problems  of  harmonic  linearization  and  frequency  methods  of 
the  investigation  of  the  SPS  are  examined  in  references  [2],  [3]  and 
[*»]. 

Let  us  examine  the  digital  tracking  system  with  variable  struc¬ 
ture  (TsSS  SPS)  in  the  control  of  pit-shaft  hoisting  machines  along 
a  track.  Let  us  take  the  closed  and  opened  structures  as  the  initial 
structures.  The  component  (Fig.  1),.  connected  to  the  input  of  the 


- 


control  system,  serves  as  the  element  for  switching  of  the  structure. 
The  logical  law  of.  switching  of  the  structure  has  the  form 


/  ',  fix)  - 


when  '/.*• 

.v  npn  x  /(  >  0; 

when  -/.v 
0  npii  x  <  o, 


(1) 


where  F(x,  px)  is  the  nonlinear  function;  x  -  the  output  value. 


Fig.  1.  Block  diagram  of  the  switching  element:  v#  ,  Vj  -  acting 
and  assigned  rates  of  movment;  hj  ,  Xj  -  traversed  and  assigned  paths 
of  movement.  Key:  1)  Rate  comparison  circuit;  2)  Coincidence  cir¬ 
cuit;  3)  Error  signal  of  track;  4)  Comparison  circuit  of  track*; 

5)  "greater";  6)  "less."  [Translator’s  note:#track  and  path  are  used 
for  the  same  word  in  the  Russian.] 

Since  the  element  for  switching  of  the  structure  is  connected  at 
the  output  of  the  digital-analog  (Ts-A)  converter  of  the  register  of 
the  matching  of  the  track  [path],  then  the  statistical  characteris¬ 
tic  and  also  the  input  and  output  signals  of  the  combined  nonlinear 
element  (element  of  switching  of  the  structure  +  converter  Ts-A), 
in  the  presence  of  a  constant  and  slowly  changing  shift  in  x°,  will 
have  the  form  depicted  on  Fig.  2.  In  this  case,  with  the  expansion 
o^  the  nonlinear  function  F(x,  px)  in  Fourier  series,  there  appears 
the  slowly  changing  constant  component,  owing  to  which  the  harmonic 
linearization  of  nonlinearity  takes  such  a  form  [5] 
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x*  is  the  self-oscillating  component  with  a  slowly  changing  amplitude 

A  and  frequency  o  ;  «  -  -  are  si  i.  ?  -  switching  angle  of  the  structure; 

-  4 

A  «  A  -  amplitude  of  the  input  signal;  S  -  quantization  in¬ 
crement  . 


Fig.  2.  Statistical  characteristic  and  nature  of  the  input  and  out¬ 
put  signals  in  the  SPS  with  the  Ts-A  converter. 

If  we  consider  that  at  the  moments  of  dimensionless  time  and 
*f«  the  output  function  F(x,  px)  by  a  jump  increases  by  one,  and 
at  moments  fit  and  iff/  by  a  jump  is  decreased  by  one,  and  at  other 
time  values  it  remains  constant  and  with  a  change  in  the  angle  from 
Pq  to  Tt  and  from  ^  to  2  K  is  equal  to  0,  then  the  coefficients  of 
harmonic  linearization  F°,  a  and  b,  computed  from  formulas  (3),  (4) 
and  (5),  have  the  form 
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where  i  is  the  integer  of  steps  of  the  characteristic  5  k  -  maximal 
number  of  steps  of  the  characteristic  covered  by  the  input  signal. 

For  the  TsSS  of  the  shaft-pit  hoist  hoisting  machine,  the  multi- 
step  characteristic  can  be  reduced  to  a  relay  by  applying  the  variable 
quantization  increment  [6].  By  keeping  this  condition  in  mind  and 
knowing  that  in  the  slip  mode  the  control  action  of  the  relay  element 
is  a  sequence  of  pulse  of  one  sign  [73,  let  us  write  conditions  of 
the  existence  of  the  slip  mode  for  one  quantization  increment. 


Since  the  condition  of  operation  in  one  direction  has  the  form 

A  f-  0.5  V-  | k°\  0.5  -  a\  (9) 

and  the  condition  of  the  absence  of  operation  in  the  other  direction 
and  in  the  direction  of  the  change  in  the  quantum  is 


Ja''|  r0.5  >  ec.£<  U°|  '  0,5, 

1,5-- U'|  >  A,  eiJni  0,5s!.vc|-  1,0, 


(10) 


then  the  condition  of  the  existence  of  the  slip  mode  denotes  the  com¬ 
bined  fulfillment  of  conditions  [9)  and  (10). 

Taking  conditions  (9)  and  (10)  into  account,  coefficients  of 
harmonic  linearization  for  the  relay  characteristic  (i  =  1)  in  the 
slip  mode  will  take  the  form 
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By  representing  the  slip  movement  as  a  slow  change  in  the  coord¬ 
inate  and  the  pulsations  imposed  on  it,  we  will  distinguish  the  slowly 
changing  and  periodic  components.  Then  instead  of  one  harmonically 
linearized  equation,  we  will  obtain  two  kinds,  respectively,  for  the 


slowly  chancing  and  self-oscillating  components  of  the  slip  movement 

/</>).<•’  i  °  S'  r  f  /j; 

/  /I-..V*  r  \\t>)  (.!+“/').«*  0.  ^  ^ 

where  L(p),  N(p)  and  S(p)  are  polynomials,  and  /(t)  is  the  external 
disturbance. 


To  find  the  periodic  solution,  let  us  examine  the  second  equation 
of  system  (1^).  The  differential  equation  of  the  TsSS  SPS  by  the 
pit-shaft  hoist  for  the  self-oscillating  component  will  have  the  form 


where  coefficients 

Ar~~<hl  A,  <Vta,s;  /'3-*0,0S<j,*:  44-0,008«V 

y, 

satisfy  the  criterion  of  optimal  figure  of  merit  [6]; 

the  time  factor;  -  calculated  duration  of  the  transient  process. 


Having  substituted  p  =  Jt>  into  equation  (15),  let  us 
the  real  and  imaginary  parts  and  equate  them  to  zero 


X\A<  jc°, 
>'{  A,  x\ 


i*,  a'i  A  ,«>*  —  A Jf®,  u>.  a)~0; 

.  «  ,  l'  .i  ,  A  «>,  •**. 

d,  a,---  to.  —  /i,ui5  4-  /l,  |-  A:, - — -  P* 


select 


(16) 


We  will  solve  the  obtained  system  of  equations  grapho-analytical- 
ly.  We  define  the  amplification  factor  from  the  the  first  equa¬ 
tion  of  system  (16) 


ai'A,  a®,  u>,  »i 


(17) 


To  find  the  amplitude  of  pulsations  A  and  frequency  o> ,  let  us 
use  the  following  method.  From  the  system  of  equations  (16),  we 


find 


A  *°.  1,1 .  O  _  A4 

a(A ,  x  ,  to,  *)  (18) 


Equation  ( 1 8 )  can  be  solved  graphically  (Fig.  3)*  To  do  this, 
from  equations  (12)  and  (13)  we  construct  in  the  second  quadrant 
dependences  of  quantity  d  on  the  amplitude  of  pulsations  A  for  dif¬ 
ferent  shifts  in  x°  when  a  =  const.  From  equation  (12)  we  plot 
the  dependence  of  the  coefficient  of  harmonic  linearization  &  on  the 
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M  a 

amplitude  of  pulsations  A  for  different  values  of  shifts  in  x°  at 
the  same  value  of  «  ,  and  from  equation  (18)  we  plot  the  cependenco 
of  the  ratio  of  d  on  the  frequency  of  pulsations  y  _  .w.  ,  Being 

assigned  values  of  <*V  ,  it  is  possible  for  each  taken  values  of  * 
and  x°  to  determine  the  amplitude  of  pulsations  A  and  coefficient 
of  harmonic  linearization  ,  and  then  from  formula  (17)  values  of 
the  necessary  amplification  factors  of  the  TsSS  SPS  in  the  slip  mode 
are  found.  The  sequence  of  the  determination  of  A  and  is  shown 
by  arrows  on  Fig.  3,  and  the  dependence  of  the  amplification  factor 
on  the  amplitude  of  pulsations,  found  by  the  grapho-analytical  method 
described  above,  is  given  on  Fig. 


Fig.  3*  Sequence  of  determination  of  (A,  x°,  c*>  ,  «  )  according  to 
the  assigned  w  ,  ft.  and  x°- 


Fig.  Dependence  of  the  amplification 
factor  of  the  TsSS  SPS  on  the  amplitude 
of  pulsations  in  the  slip  mode:  -  -  -  - 
-  unstable  part  of  the  characteristic; 
-  -  stable  part  of  the  charac¬ 
teristic  . 
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Fig.  5.  Smoothed  functions  of  the  shift 
in  the  nonlinear  element  for  the  slowly 
changing  component  in  the  slip  mode. 


In  order  to  determine  the  slowly  changing  component  F°(x°),  in 
equation  (11)  it  is  necessary  to  exclude  A,  for  which  we  plot  the  de¬ 
pendence  of  the  amplitude  of  pulsations  A  on  the  shift  in  x°  when 
Ajj  s  const,  having  used  the  graph  depicted  on  Fig.  4,  taking  condi¬ 
tions  (9)  and  (10)  into  account.  Then,  according  to  formula  (11), 
exluding  A,  we  plot  the  dependence  F°(x°).  The  found  functions  of 
the  shift  in  the  nonlinear  element  in  the  slip  mode  (when  s  0.5-rr) 
are  given  on  Fig.  5.  After  the  linearization,  the  obtained  charac¬ 
teristics  arc  described  by  equations  (when  =  0.5fO 

/  M  0,1V)  npit  A,  C,u  mq\ 

/  '  <M.;s  upa  A:,  on"  !<."».»,*: 

/•’-  O.S'.i.V  lien  A:.  I » J  * 

1 ''  II  |  Ml  A.  -Q.tHlO.ViW,' . 

Figure  6  gives  dependences  of  the  amplification  factor  -- ,  , 

calculated  by  the  described  method,  on  the  curvature  of  the  linearized 
characteristic  of  the  nonlinear  element  k  for  the  slowly  changing 
component  for  different  switching  angles  of  the  structure  cC  . 

By  using  the  obtained  dependences  (19)  and  substituting  therm 
into  the  first  equation  of  system  (14),  it  is  possible  to  construct 
by  the  method  of  trapezoidal  frequency  characteristics  the  transient 
process  for  the  slowly  changing  component  in  the  slip  mode  both  for 
the  controlling  and  disturbing  actions.  The  curve  of  the  transient 
process  for  the  slowly  changing  component  for  the  controlling  action 
when  =  0.2 rr  and  =  0,001  plotted  by  the  method  of  trupe- 

soidal  frequency  characteristics,  is  given  on  Fig.  7.  Curves  of  the 
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transient  process  for  different  values  of  the  amplification  factor 
A,  can  be  similarly  plotted. 


Fig*  6.  Dependence  of  the  amplification 
factor  of  the  system  on  the  curvature  of 
the  linearized  characteristic  of  the 
nonlinear  element. 


Fig.  7.  Transient  process  of  the  TsSS 
SPS  according  to  the  controlling  action. 

Here  ,  where  t  is  the  time. 
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Results  of  theoretical  studies  of  slip  modes  of  the  TsSS  SPS 
were  conducted  on  an  electrodynamic  model. 

Figure  8  gives  an  oscillogram  of  the  section  of  the  change  in 
rate  from  V  =  10  m/s  to  V  =  0  with  constant  acceleration  for  selected 
values  of  parameters  Ct  =  0.2  7 T,  =  O.OOla^  and  =  42.14  s”1. 

The  calculation  of  parameters  of  the  model  is  carried  out  according 
to  the  method  described  in  [6].  Operation  of  the  TsSS  SPS  was  con¬ 
ducted  in  the  slip  mode  with  a  variable  quantization  increment  .  The 
criterion  of  the  coincidence  of  results  of  theoretical  and  experi¬ 
mental  studies  is  the  magnitude  of  error  [mismatch]  of  the  path, 
which,  according  to  the  calculations,  must  not  exceed  the  quantiza¬ 
tion  increment  for  any  value  of  the  rate.  The  slip  mode,  which 

8 


consists  in  oscillations  in  the  error  of  the  path  by  a  magnitude  of 
one  quantum  is  clearly  marked  on  the  oscillogram  (Fig.  8).  Owing  to 
the  fact  that  the  traveling  quantum  is  selected  as  variable  and  in¬ 
creasing  with  an  increase  in  the  rate,  the  frequency  of  the  oscilla¬ 
tions  is  changed  insignificantly  with  a  change  in  the  rate  from  max¬ 
imal  to  zero,  Evident  are  pulsations  of  the  slowly  ohanging  com¬ 
ponent,  which  consist  in  pulsations  of  rate.  The  error  with  respect 
to  rate  is  such  that  the  magnitude  of  the  error  of  the  path  does 
not  exceed  the  quantization  increment,  which  indicates  the  high-speed 
operation  of  the  system  of  automatic  control. 


t  . 

I . -  •  •  ' 


Fig.  8.  Oscillogram  of  the  section  of  slowing  down  of  the  TsSS  SPS 
in  the  slip  mode:  1  -  assigned  rate  of  movement;  2  -  actual  rate  of 
movement;  3  -  current;  4  -  error  of  the  path. 


The  proposed  method  makes  it  possible,  with  an  accuracy  suffi¬ 
cient  for  practical  calculations,  to  calculate  parameters  of  the  slip 
mode  in  the  digital  tracking  control  system  with  variable  structure, 
the  dynamics  of  which  is  described  by  the  nonlinear  differential 
equation  of  a  high  order. 
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